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Abstract. 

Using the fluid representation, we formulate the conditions for the appearance of all four types finite-time future singularity 
in modified gravity in accelerating FRW universe. It stressed that number of standard quintessence/phantom DE theories 
(including scalar, fluid, DBI ones, etc) brings the accelerating cosmology to future singularity precisely in the same way 
as singular modified gravity DE. The viable non-singular modified gravity unifying the early-time inflation with late-time 
acceleration is considered. It is shown that adding such non-singular theory to another realistic singular modified gravity which 
has the accelerating solution with future singularity may cure the singularity of resulting combined model. This universal 
scenario may be naturally applied to standard singular DE models as well as to inflationary theories with initial singularity. 
This suggests the additional fundamental reasoning for modification of General Relativity. 
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INTRODUCTION 

Modified gravity suggests the gravitational alternative for the unified description of the early-time inflation with late- 
time acceleration and dark matter. It is remarkable that in this scenario for the universe evolution there is no need in 
the introduction of extra dark components like inflaton, quintessence, dark fluid, dark matter particles, etc. The only 
early-time and late-time gravitational action is changed if compare with General Relativity in such a way that local 
tests and cosmological bounds are not violated. Hence, the gravitational field is responsible for dark behavior of the 
universe. The comparison of different modified gravities aimed to the accelerating universe description may be found 
inref.H. 

Imagine that we consider the (accelerating) FRW universe. In this case, whatever is modified gravity, it may be 
represented as General Relativity with some (generalized) fluid. Different properties of modified gravity may be 
clarified in this representation. For instance, for realistic theories which unify the early-time inflation with late-time 
acceleration and which we consider here it is clear that they behave as effective phantom/quintessence fluids at the late 
times. Then, such modified gravities may bring the evolution to the finite-time future singularity of one of four types 
according to the classification of ref.J2l in the same way as other phantom/quintessence dark energy models. There is 
no any qualitative difference between modified gravity and other dark energy models in this respect. (Reversing the 
time flow, the finite-time future singularity is nothing else as initial singularity for number of inflationary theories.) 

In the present paper we explicitly demonstrate that the singularity appearance is typical property for number of 
modified gravity with the effective quintessence/phantom late-time behavior precisely in the same way as for other 
(scalar/fluid/composite) dark energy models. The appearance of all four types of future singularity in unified modified 
gravity occurs as it was demonstrated in ref.yfl (for first observation on Big Rip singularity in modified gravity and its 
avoidance by adding of /? 2 -term see ref.[4, 5]). The realistic non-singular modified gravity which unifies the early-time 
inflation with late-time acceleration may be considered as the tool to prevent the future singularity. In particular, it is 
demonstrated that adding such non-singular theory to another realistic modified gravity which has the accelerating 
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solution with future singularity may cure the singularity appearance in combined realistic modified gravity. In the 
same way, the initial singularity may be cured unless the classical description is valid. 



THE FORMULATION OF MODIFIED GRAVITY AS GENERAL RELATIVITY WITH 
GENERALIZED FLUID AND FINITE-TIME FUTURE SINGULARITIES 

Let us start from the general modified gravity with the action: 

S = jd 4 x^^^^R + f^R,R^\R^ va pR^ a P,aR,a- l R,---^+L m Y (1) 

where all combinations of local and non-local terms are possible, L m is matter Lagrangian and the function f(R, ■ ■ ■) 
may also contain gravitational partners (say, dilaton, axion, etc. in string-inspired gravity). 

We consider spatially-flat FRW universe with scale factor a(t). In all cases for theory (|T), it is possible to write the 
gravitational field equations in the form of standard FRW equations with effective energy-density p e ff and pressure 
Peff produced by the extra gravitational terms F(R,- ■ • ) and L m . 

For instance, when / = f(R), one gets 

Peff = ^^f(R) + 3{H 2 +H)f\R)-lS(AH 2 H + HH)f(R^j+ Pmatter , (2) 
Peff - ^ (if ( R ) ~ ( 3//2 + H)f(R) +6 (M 2 H + 4H 2 + 6HH + H) f" (R) + 36 (4HH + ti) 2 f"'(R)^j 

+Pmatter ■ (3) 

For details of explicit fluid presentation of f(R) gravity, see |@]. In case of Gauss-Bonnet modified gravity f(G) 10]: 
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In the same way one can get the effective gravitational pressure and energy density so that the equations of motion for 
arbitrary modified gravity can be rewritten in the universal FRW form typical for General Relativity: 

^// 2 =Peff, Peff = -^T (2// + 3/Z 2 ) . (5) 

There are just standard FRW gravitational equations. 

Formally, the modified gravity equations © have just the same form as for GR with matter. The only difference is 
that the effective energy-density and pressure are caused by extra gravitational terms due to the modification of the GR 
Lagrangian. It is quite well-known that equations of motion © can describe the accelerating (early-time or late-time) 
epoch when all or part of the (effective) energy conditions are violated. For instance, when w e s < — 1, one gets the 
effective phantom (super)acceleration evolutionary phase (all energy conditions are violated). When — 1 < w e ff < — 1/3 
(Peff = WeffPeffX one gets the effective quintessence acceleration epoch, in this case, only part of energy conditions is 
violated. For w s s = — 1, the de Sitter acceleration evolution emerges. 

The remarkable property of phantom/quintessence accelerating universe is the occurrence of the finite-time future 
singularity which is caused by the violation of the energy conditions. The phantom evolution when it is not transient 
always ends up with so-called Big Rip singularity. The quintessence accelerating evolution may enter (or may not 
enter) to soft finite-time future singularities. The classification of four possible finite-time future singularities for late- 
time accelerating universe has being made in ref . 12|] . 

The clarifying remark is in order. Big number of phantom/quintessence dark energy models (say, fluid, scalar, spinor, 
DBI, tachyon, etc.) brings the evolution to the finite-time future singularity. From the above, equivalent description 
it is clear that number of modified gravities also describes the effective phantom/quintessence accelerating universe 



which evolves to finite-time future singularity. It is important to understand that the singularity occurrence is general 
feature of the phantom/quintessence dark energy models which may be fluid, scalar, or modified gravity. Moreover, 
even the classification of finite-time singularities for modified gravity is given by the same scheme [2] where p and p 
should be substituted by p^, p e ff : 

• Type I ("Big Rip") : For t — > t s , a — ► oo, p eff — > oo and \p e g\ — ► °°. This also includes the case of p e ff, p e s being 
finite at t s . 

• Type II ("sudden") : For t —> t s , a — ► a s , p e ff — > p s and |/? e ff | —* °° 

• Type III : For t — > f s , a — > a s , p e ff — > oo and |/> e g | — > oo 

• Type IV : For f — > f s , a — * a i; p e ff — > 0, |/> e ff| — > and higher derivatives of H diverge. This also includes the case 
in which p e g (p e ff) or both of /? e ff an d Peff tend to some finite values, while higher derivatives of H diverge. 

Indeed, it has been observed the approach to Big Rip singularity in the modified gravity in ref.[4]. However, it was 
demonstrated there that phantom phase becomes transient due to presence of 7? 2 -term and no singularity occurs. In 
ref.Jll, it was clearly demonstrated the appearance of Big Rip type singularity in several f(R) modified gravity models 
(see also examples of modified gravity which describe the phantom era and occurrence of Big Rip singularity in |Q]]). 

The explicit construction of f(R) gravity models which lead to all four possible types of future singularity after 
the corresponding quintessence/phantom era has been presented in ref.yfl. Even more explicit models of f{R) and 
other modified gravities which realize the four possible finite-time future singularity have been presented in ref.||8|. 
The search of singular modified gravity models was done using the reconstruction scheme developed in refs.Jgt. In 
this way, the modified gravity which realizes finite-time Type I,II,III or IV future singularity may be constructed in 
the same way as for traditional DE models (scalar, vector, spinor, etc). Moreover, it was shown in refs.JH |H [Toll , 
how to make the non-singular model from the initially singular modified gravity. Explicitly it was considered the 
addition of 7? 2 -term (as in the non-singular modified gravity model 10]), or some other gravitational terms |§] relevant 
at the early universe, or the account of the quantum contribution due to conformal anomaly as scenario for curing 
the future singularity of any type. (For related discussion of Type II finite-time future singularity in specific model 
of f(R) gravity and its avoidance by the same trick 10] of adding R 2 -term, see refs. lfTTll ). It is interesting to note that 
above four types of future singularities in FRW space-time reappear in spherically symmetric space-time as energy 
density and pressure singularities at finite radius IU2I1 . It is precisely these singularities which may lead to problems 
in the formation of relativistic stars and black holes in any dark energy model evolving to future singularity. In the 
next section, we discuss the realistic f(R) gravity which is non-singular and which unifies the early-time inflation with 
late-time acceleration. 



NON-SINGULAR REALISTIC f(R) MODEL UNIFYING THE EARLY-TIME 
INFLATION WITH LATE-TIME ACCELERATION 



, . aR + - BR 
f(R) = t — • (6) 



Let us consider the following f(R) model 1FJ 



For the case m = l = n, the model has been analyzed in ref.[ 13]. It was shown that the model may describe the realistic 
unification of early-time inflation with late-time acceleration (for a first proposal of the modified gravity with such 
unification, see fl4ll ). For the realistic unification, By/ a ~ 10 228 (" _1 )eV 2 in the solar system. This gives extremely 
small (non-observable) correction to the Newton law if n > 2. In this case, there is no the observable fifth force (the 
corresponding force is much below the observational bounds). It is not difficult to check that for « = 2, the above 
model © which passes the local tests is also non-singular. Hence, the model © suggests the generalization of the 
scenario to cure the future finite-time singularity via adding the /? 2 -term [4]. Actually, the same mechanism is used: 
one can consider the addition of the term R 2 f(R) where lim^o/CK) = c h li m «^~/(^) = c 2 as the scenario to remove 
the future singularity. In principle, other higher derivative terms may be also proposed for this purpose. Hence, one 
can consider the realistic dark energy alternative f(R) gravity from the class suggested in ref . II 1 511 (or more general 
realistic models [ 16] unifying the early-time inflation with late-time acceleration). 

In case that such theory contains the finite-time future singularity (as it happens in many other dark energy models 
|0]), one can add /? 2 -term [4] or R 2 f(R) (see above) in order to cure the finite-time singularity. Of course, one can also 



take into account quantum effects due to conformal anomaly 10, [8| or quantum gravity effects lfl7ll or generalized fluid 
of special form HE] to prevent the occurrence of future singularity. Moreover, the use of the function of the form 
© to cure the finite-time singularity seems to be the most promising one. The reason is that the terms of such sort are 
mainly relevant at the early-time era, so eventually they do not spoil the local tests of the theory. From another side, 
such terms give the contribution to the cosmological perturbations so probably they may be used to make that more 
consistent with the observational data. 

As the realistic example, we consider the viable theory [16]: 

(R-R ) 2k+l +R 2k+l 
f(R) = y — f ^ 5 (7) 

/o + / 1 {(/?-/?o) 2i+1 +< +1 } 

It has been shown lfl6Tl that for k > 10 such modified gravity passes the local tests. It also unifies the early-time inflation 
with dark energy epoch. 

In (Q, i?o is current curvature Rq ~ (lO~ 33 eV) 2 . We also require, 

Rl n 1 

f °~ir ■ /i = at- (8) 

Here A, is the effective cosmological constant in the inflation epoch. When R 3> A,-, f(R) Q behaves as 

1 , /o 



We now use the trace equation, which is the trace part of the FRW gravitational equation 

3Df (R)=R + 2f(R) -Rf'(R) - k 2 T . (10) 

Here T is the trace of the matter energy-momentum tensor. Then assuming the FRW metric with flat spatial part, one 
finds 

fl~(? -?r 2/(2 " +3) , (ii) 

which diverges at finite future time t = to- By a similar analysis, we can show that if f(R) behaves as f(R) ~ R a for 
large R with a constant a, a future singularity appears if a > 2 or a < 0. Conversely if 2 > a > 0, the singularity does 
not appear. Then by adding the previous term R 2 f(R), where lim^o/(^) = c i, ^ m R^oof(R) = C2, to f(R) in the 
future singularity ( fTTT i disappears. 

Let us observe the above situation in more detail. Now we assume 

f(R)^F + F l R a , (12) 

when R is large. Here Fq and Fx are constants where Fq may vanish but we assume F\ ^ 0. In case of (O, we have 

^o = -^, *i = 4, a = -(2« + l). (13) 
/i ft 

Under the assumption (fT2l . the trace equation (TTOb gives 

3F 1 DR"- 1 = l% , rpa ^en a <0 or a = 2 

| (2 - a)FiR a when a > 1 or a ^ 2 v ' 

In the FRW background with flat spatial part, when the Hubble rate has a singularity as 

ho 

with constants ho and /3, the scalar curvature R = 6H + 12H 2 behaves as 

^ when^>l 
R~l ^±i^ whenj3 = l . (16) 
(^ when^<l 



In <TT3T > or ( TToT ). j3 > 1 case corresponds to Type I (Big Rip) singularity, 1 > /3 > to Type III, > /3 > — 1 to Type II, 
and j3 < — 1 but j3 7^ integer to Type IV. By substituting ([Tol l into (fT4l i. we find that there are two classes of consistent 
solutions. The first solution is specified by j3 = 1 and a > 1 but a^2 case, which corresponds to the Big Rip 
(ho > and t < to) or Big Bang (ho < and t > to) singularity at t = to. Another one is a < 1, and j3 = —a/ (a — 2) 
(— 1 < j3 < 1) case, which corresponds to ( fTTT i and to the II Type future singularity. In fact, we find a = —2n — 1 and 
therefore J3 + 1 = —2 / (2n + 3). We should note that when a = 2, that is, f(R) ~ /? 2 , there is no any singular solution. 
Therefore if we add the above term R 2 f(R), where lim^o/C^) = c\, lim^oo/(,R) = c%, to f(R) in (0, the added term 
dominates and modified f(R) behaves as f(R) ~ R 2 , the future singularity ( fTTT i disappears. We also note that if we add 
R"-term n — 3,4,5, • • • , the singularity becomes (in some sense) worse since this case corresponds to a = n > 1, that 
is Big Rip case. Using the potential which appears when we transform F(/?)-gravity to scalar-tensor theory [3], it has 
been found that the future singularity may not appear in case < a < 2. 

Thus, we demonstrated that combination of two realistic modified gravities where one of them is non-singular one 
may heal the finite-time future singularity of combined realistic model. On the same time, the nice properties of the 
models: the successful passing of local tests as well as unified description of the inflation with dark energy epoch 
remain to be the same. 



DISCUSSION 

In summary, using the fluid representation, we demonstrated that future singularity appearance is typical property 
of number of modified gravities with the effective quintessence/phantom late-time behavior in the same way as for 
other simpler (scalar/fluid) dark energy models. These four different future cosmological singularities may manifest 
themselves as effective pressure/energy-density divergence at finite radius in spherically-symmetric space. This may 
cause number of problems with black holes and relativistic stars formation for any singular dark energy (as is shown 
in IU2I1 ). not only for alternative gravity dark energy models. Moreover, the easiest way to cure the future singularity is 
again to call for the modification of gravity, for instance, adding the /? 2 -term which is relevant at the very early universe. 
Hence, starting from singular fluid/scalar dark energy one is forced to modify its gravitational sector in order to cure 
the singularity. Is not then more natural to start from the modified gravity from the very beginning? Moreover, the 
combination of non-singular realistic modified gravity with singular realistic model unifying the early-time inflation 
with late-time acceleration leads to the resulting realistic non-singular theory. We concentrated mainly on f(R) gravity 
but other modified gravities may show very similar singular behavior for effective quintessence/phantom late-time era 
as is explained in ref.[8]. The same universal tools to cure the future singularity may be used for any modified gravity. 
Finally, it is interesting to mention that above scheme easily accounts for description of DM effects within the same 
modified gravity (for review of DM properties from modified gravity, see I19ID . 
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